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Abstract 

A theoretical analysis of surface fluctuations of a freely standing thermotropic 
smectic-A liquid crystal him is provided, including the effects of viscous hydro- 
dynamics. We hnd two surface dynamic modes (undulation and peristaltic). 
For long wavelengths and small frequencies in a thin him, the undulation 
mode is the dominant mode. Permeation enters the theory only through the 
boundary conditions. The resulting power spectrum is compared with exist- 
ing experiments. It is also shown that feasible light scattering experiments on 
a freely standing smectic-A him can reveal viscosity and elastic coefficients 
through the structure of the power spectrum of both the undulation and peri- 
staltic modes. 
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Recent experiments on dynamics and instabilities of soap |IJ and smectic-A films have 
raised interest in the general question of dynamics of freely-standing films. The interactions 
between bulk elasticity and surface tension make freely-standing smectic-A films (FSSF) 
suitable systems for studying finite-size and surface effects; hence, fluctuations in FSSF are 
an important subject for both theoretical and experimental study [^]. However, during the 
past decade only static properties of these systems have been considered. To our knowledge, 
the only experimental work on dynamic light scattering of FSSF was carried out without 
systematic theoretical analysis. Hence a theoretical investigation on the dynamics of smectic- 
A liquid crystals in the presence of free surfaces can help us to gain deeper insight to the 
physics of FSSF and to provide a basis for understanding new experiments. 

At free surfaces, smectic layers always orient parallel to the smectic/air interface ||. 
When the system is driven out of equilibrium, surface and bulk elasticity and hydrodynamic 
effects give the film very complicated dynamical properties. In this paper we give the first 
theoretical calculation of the power spectrum of the surface fluctuations by studying the 
linearized hydrodynamic equations developed by Martin et al. |J. We show that for thin 
films, the long wavelength surface fluctuations have two modes. The undulation mode is the 
dominant mode at low frequency and long wavelength. The peristaltic mode is expected to 
be very small. Only when the thickness of the film becomes very large does the peristaltic 
mode become comparable to the undulation mode. Permeation processes are important 
within the boundary layer |7|,|8| in order that the proper boundary conditions are satisfied. 
However, our calculation shows that the power spectrum of the surface fluctuations are 
independent of the permeation constant. We also show that the earlier experimental study 
of dynamic light scattering of FSSF |4j] was performed in the limit of a very thin film such 
that the power spectrum has the same form as that of a soap film, i.e., bulk elasticity makes 
no contribution to the power spectrum. The scaling relations suggested in Ref. [Q for the 
peak position are not valid in general, especially when the thickness of the film is increased, 
thus allowing the bulk elasticity of the smectic material to play a role in the dynamics. For a 
reasonably thick film, the power spectrum can have, in addition to a single undulation peak, 
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some additional structure which reveals the interaction between the two free surfaces and 
the contribution from bulk properties. We suggest that more light scattering experiments 
on a FSSF will be great help for an understanding of the interplay of the surface and bulk 
properties in this layered system. 

In the ground state, an ideal smectic-A phase consists in a uniform piling of planar, 
parallel and equidistant layers of molecular thickness. We use a continuum description and 
take the average layer normals parallel to the z-axis. When small fluctuations are present, 
the bulk elastic free energy of smectic-A phase is given by @ 



1 I „ / du\ 2 „ / d 2 u d 2 u K 2 



where u(r, t) is the layer displacement from the equilibrium position at position r at time t, 
B and K\ are, respectively, the the layer compression and undulation elastic moduli. The 



characteristic length A = yK\/B is typically of the order of the layer spacing (~ 10 7 cm) 
0- 

The equations for viscous flow in smectic A liquid crystals were first written down by 
Martin et al @. In the absence of topological defects, the system satisfies the equations of 
motion in bulk 

Ov' 

p-r^- = -dip + djcr'ij + h 5 iz (2) 



and 



du . . 

= v z + C P h , (3) 

where a' is the viscous stress tensor, ( p is the permeation constant and h is defined by 

In (0) we sum on repeated indices, and dj = d/dxj. An important length scale associated 
with permeation is the boundary layer 5. Within this distance to the boundary, permeation 
takes place to satisfy proper boundary conditions of the system under consideration |7| -|9|. 



We will show that the boundary layer is associated with the "permeation modes". The 
boundary conditions for free surfaces are discussed below. 

We consider a freely standing Smectic-A liquid crystal film with film normal in the in- 
direction. As noted, in equilibrium, the layers are parallel to the free surfaces. When the 
surfaces are perturbed by an external force, elastic forces and dissipative effects act to drive 
the system back to the equilibrium configuration. In this letter we study such a system by 
including the hydrodynamics in the film. 

The geometry of the film is shown in Fig. 1; the film extends from z = —d to z = 
but is otherwise without boundaries in the x, y directions. We consider surface light 
scattering experiments with momentum transfer vector q in the ^-direction; we further 
assume translational invariance in the ^-direction. The displacement of the upper and lower 
free surfaces from their equilibrium values are described by two functions ( + {x) and C ^)- 

We look for surface wave solutions with the form 

v z = { A t e Skqz + A^e- Skg(z - d) } e lwt+iqx (5) 

k 

where Re(Sfc) > 0. The + (— ) modes are the upper (lower) surface modes which have their 
maximum amplitudes at z = (z = —d). 

Two combinations of viscosities rf and 773 || and a length scale /t _1 = ^CpVs ~ 10~ 7 cm 
|7]-|| enter the analysis. A dimensionless frequency is conveniently defined as Q = The 
characteristic frequency Kiq 2 /i] 3 for the decay rate of the bulk undulation mode, reduced by 
r/ 3 q 2 /p, yields a dimensionless frequency fx = We assume incompressibility (the motion 
of the fluid is slow compared to the propagation of sound), long wavelength ,i.e., (n/q) 2 ^> 
1, Xq I. We consider low frequency satisfying conditions ^(Ag) 2 1, y-(Ag) 2 I. 
This means we consider a regime in which the velocity of the surface wave is slow compared 
to the typical velocity of "second sound", \J B / p [0. Furthermore, since then the frequency 
cannot be as large as Ki\~ 2 /rj3, permeation cannot have a significant contribution to the 
bulk undulation mode. J7| In this frequency and wavenumber domain we find three spatially 
decaying modes in the z-direction. One of them has long spatial relaxation length and 
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satisfies 



Si 



(Ag) 2 



l - £(i 



n) 



1 + - (Ag) 2 (fi-^-2 

M V »?3 



= (Ag) 2 [/(^,g)f 



(6) 
(7) 



(8) 



The other two modes (Si and S2) are large compared to unity and satisfy 

+ 9. ( K \fs 2 + (-X = 0. 

A* W 

In typical light scattering experiments g ~ 10 2 — 10 4 cm -1 , to < 10 8 rad/s. For typical 
materials % ~ 1 p so that the long wavelength, low frequency conditions above are satis- 
fied. Notice that when the denominator of the function /(O, q) defined above is significantly 
different from unity, we cannot separate the solutions for Si, S2, S3 in the way mentioned 
above. The mode with relaxation S3 is related to the dynamic generalization of the pene- 
tration length of layer distortion in the bulk due to surface undulation. [|7|pTJ| Following the 
terminology of Rapini, JTJJ from now on we refer to the S^-mode as the elastic mode. Since 
Si, and S2 depend on the permeation constant, ( p , they are called permeation modes. For 
e Siqd , e S2qd 3> 1, i.e., d^> (Siq)^ 1 , (S 2 q)~ 1 ~ 5, the film thickness is large compared to the 
boundary layer 6 (identified as the exponential decay length), and the permeation modes 
relax in the bulk. 

The surface displacements ( + , (~ satisfy the boundary condition for the normal compo- 
nent of the velocity in linear theory 



dt 



V 



z\z=0(-d)- 



(9) 



The other boundary conditions can be understood from the covariant elasticity theory of 



smectic-A developed by Kleman and Perodi [11|. For free surfaces the normal components 



of the stress tensor as well as the normal component of the permeation force |L2[ should 
vanish. The condition that a xz = on the free surfaces reveals A% ' <C ^.3 f° r k — 1, 2, 
which indicates that the permeation modes have negligible contributions to the dynamics 
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of the system. The condition that a zz = yields the following relation between the surface 
displacements and the external forces P+ t , P~ xt which are assumed to act on the free surfaces, 



P +(-) 

r ext 



du <9 2 c +H 

P- (X zz + B^r + a- 



(10) 

z=0(-d) 



dz dx 2 

where a is the air-film surface tension. The normal component of the permeation force, in 
the system under consideration, is 

at z — 0, —d for all x. The linear response function X connects the surface displacements 
with the external forces through 

C(q,w) = -X(q,w)P ext (q,w)A , (12) 

where A is the surface area, P ex t are vectors for surface displacement and external force 
respectively, e.g., X is a 2 x 2 matrix for the response function. 

Putting Eqns. (|TT|) , (|I0|) , and fll2"|) together our calculation leads to two surface dynamic 
modes; they are, respectively, the undulation mode with amplitude 

c u = \{c + + C) (13) 

and the peristaltic mode with amplitude 

C p = \ (C + - C") • (14) 

The response functions are in turn given by 

J£U _ 1 1 

2a? A l + ^/(n, g )tanh(^/(n,g)) ' 

and 

X P = 7^r : . - ; \ T^— , (16) 



2aq 2 A 1 + g f(Q, q ) co th(^f(Q,q 
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where g = \JBK\j a measures the relative importance of bulk elasticity and surface tension. 
When the contribution from bulk elasticity and viscosity dominates, terms involving g are 
dominant. These equations are the central result of our calculation. 

In the limit of infinite thickness, d — > oo, the response functions of both modes become 
the same, indicating that the interactions between the two surfaces vanishes for large thick- 
ness. For low frequency (Q <C 100) and infinite thickness, the result agrees with the earlier 
calculation by Rapini [TO] for an infinitely thick smectic-A material with a free surface. 



We now introduce a natural frequency of surface motion as Uq = ^q 2 , and dissipation 
coefficient 7 = In the range of the only available experiment @ (thin film, d ~ 100 

nm, small wavenumber, q ~ 10 4 cm -1 ), Xq 2 d f(Q,q)/2 <C 1. The response function is 



approximately 



1 + lg [mqWXqZd 

~ {ul-u^+iw ' (17) 
which is the same as the response function of a damped driven simple harmonic oscillator. 
The peak position in the range of weak damping {^-oj <C \f^i ( l ) * s 

1 2a . . 

This is the same as for a soap film and gives the scaling relation which the experimental data 
satisfy. In general we may not have a system with weak damping, and the peak position 
should be slightly modified. On the other hand, the width of the undulation mode peak 
iV3q 2 /p) does provide information about the viscosity 773. 

The peristaltic mode changes the layer spacing much more significantly than the undu- 
lation mode; hence we expect that the peak in the peristaltic power spectrum occurs when 
the term proportional to g in Eq. (jlBD dominates. Also in the range that we are interested 
in, the relation Q/fi ^> 1 is valid. A straightforward calculation leads to an estimate of the 
peak position of the peristaltic mode : 

Xq 2 d A(fi,g)«7r (19) 



where fi(Q,q) is the imaginary part of the function / we have defined via Eq. (|7|). In this 
approximation the surface tension plays no role in peristaltic mode, the peak position is 
essentially independent of g and has the form 

O • 2/17T 2 

(20) 

For typical material parameters the above approximation yields the peristaltic peak within 
5%. From this relation we can estimate the bulk elastic coefficients K\ and B by fitting the 
value /i and A from the peristaltic peak. 

One may ask whether it is possible to observe the peristaltic mode and the special features 
of a smectic-A liquid crystal in a free-standing film experiment. Figure 2 shows the power 
spectrum of one free surface ,i.e, ( u + ( p , for a typical choice of parameters with increasing 
film thickness. As the film gets thicker we observe that a peak develops in the higher 
frequency range. This peak is due to the peristaltic mode. When L = qd > 12, % Q > 3 
there is some extra structure arising from sources other than the peristaltic mode. Figure 3 
shows the power spectrum for the same choice of parameters with qd = 15. The contribution 
from both undulation and peristaltic modes are plotted in long and short dashed lines, 
respectively. We find that the power spectrum of the peristaltic mode is about one order 
of magnitude smaller than the contribution of the undulation mode. Notice that there is 
a second peak of the undulation mode in the high frequency region. It comes from the 
oscillating behavior of the hyperbolic tangent function when the argument ( Xq 2 d f(fl, q)/2 
) is complex. This extra structure in the power spectrum at high Q is a result of the 
interplay of bulk elasticity and the existence of the two free surfaces; it is a special feature 
of a structured fluid. Hence we conclude that for reasonable choices of material parameters 
and for an experiment with typical dynamic range, it should be possible to observe the 
peristaltic mode. Combinations of material parameters can be extrapolated from the shape 
of the undulation peak and the peristaltic peak. However, the detailed shape of the power 
spectrum is sensitive to the specific material parameters used in a laboratory experiment. 

In conclusion, we have derived the power spectrum of a freely standing smectic-A film 
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within linear hydrodynamics and assuming the absence of topological defects. The dynamics 
of this system are dominated by the elastic mode and the permeation constant does not show 
up in the power spectrum of the surfaces. The permeation process is important near the 
surfaces and helps the system to satisfy proper boundary conditions. When the thickness 
of the film is small enough, bulk elasticity does not contribute to the undulation mode, 
and the peristaltic mode is not observable. However, for a reasonably thick film the power 
spectrum does show the interplay of surface tension and bulk elasticity. The extra structure 
in the power spectrum is a special feature due to the existence of two free surfaces and the 
contribution from the bulk elasticity. We suggest that further experimental work on FSSF 
over a wider range of film thickness and wavenumber can observe these interesting features. 

We thank Prof. X-l. Wu and D. Dash for very helpful discussions, and Prof. Takao Ohta 
for his interest and assistance. H.Y.C. acknowledges financial support from University of 
Pittsburgh as an Andrew Mellon Predoctoral Fellow. D.J. is grateful for the support of the 
NSF under DMR9297135. 
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Figure Caption 

Figure 1. Schematic of a freely standing smectic-A film of thickness d. 

Figure 2. The power spectrum (natural logarithm) for a typical choice of parameters 
with increasing thickness of the film: \q = 1CT 3 , /i = 10~ 4 , g — 0.1, rj' — 773. L = qd is 
dimensionless. As the thickness increases, we can easily see the peristaltic mode and extra 
structure of the undulation mode provide additional contribution to the power spectrum. 

Figure 3. The power spectrum (natural logarithm) for a typical choice of parameters: 
\q = 1CT 3 , qd = 15.0, /i = 10~ 4 , g — 0.1, rj' — 773. The long dashed line is the undulation 
mode, short dashed line is the peristaltic mode, solid line is the total power spectrum. Notice 
that there is a small peak in the undulation mode at i Q ~ 4. 
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